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the  dislocation  structure  in  Ir  and  Rh-based  Ll2  alloys.  We  determined  the  temperature  dependence  of  the  lattice  misfit  parameter  in  y/y' 
superalloys,  and  established  its  simple  relation  with  the  alloy  phase  diagram.  We  studied  the  structure  of  dislocations  in  bcc  Nb,  Ta,  Mo  and 
W  alloys,  analyzed  the  effect  of  alloying  with  3 d  and  5 d  transition  metal  elements  on  dislocation  energetics,  and  established  the  electronic 
origins  of  the  solid  solution  softening/hardening  phenomena.  We  investigated  the  effect  of  topologically  close  packed  phase  formation  on 
the  solubility  of  interstitial  impurities  in  Cr-  and  Mo-based  alloys  and  demonstrated  their  effect  in  preventing  the  formation  of  embrittling 
phases.  We  showed  that  Co-  and  Y-based  B2  intermetallics  are  intrinsically  ductile,  established  the  dislocation  structure  in  Co-based  alloys, 
analyzed  the  origins  of  their  unusual  mechanical  behavior,  and  identified  the  microscopic  origins  of  the  yield  stress  anomaly. 
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Objectives 


The  objective  of  this  research  is  to  investigate  the  microscopic  mechanisms  governing  the 
deformation  and  fracture  behavior  of  ordered  intermetallic  aerospace  alloys  and  metals  in 
order  to  contribute  to  the  development  of  a  fundamental  basis  for  computer-aided  alloys 
design.  The  most  important  and  challenging  component  of  our  research  is  to  bridge  the  gap 
between  a  microscopic  quantum-mechanical  description  of  the  chemical  bonding  and  the 
mesoscopic  phenomena  which  govern  the  mechanical  response  of  intermetallics.  The 
emphasis  on  the  computational/simulation  approach  focuses  on  understanding  “real'’ 
materials,  which  have  an  abundance  of  “defects”  including  impurities,  dislocations  or  other 
faults,  second  phase  precipitates,  etc.  -  all  of  which  are  governed  on  the  microscopic  level  by 
the  electronic  structure. 


Research  findings 


7.  Fundamental  parameters  of  brittle-ductile  behavior  of platinum  group  metals  alloys  and 
prediction  of  cleavage-like  crack  propagation. 

Ir-  and  Rh-based  alloys  attract  considerable  attention  as  promising  materials  for  ultra-high- 
temperature  applications.  However,  polycrystalline  Ir,  Rh  and  their  alloys  normally  exhibit 
intergranular  fracture  with  limited  ductility  in  tensile  tests  up  to  elevated  temperatures  [I], 
The  lack  of  understanding  of  the  fundamental  characteristics  which  control  brittle-ductile 
behavior  in  these  alloys  has  hindered  progress  in  this  field. 

The  brittleness  of  Ir  and  Rh  is  intrinsic  and  is  not  caused  by  impurities  at  grain  boundaries. 
As  we  showed  in  [2]  the  reasons  for  this  behavior,  which  distinguishes  Ir  and  Rh  among 
otherycc  metals,  is  connected  with  a  “pseudo-covalent”  character  of  the  chemical  bonding  in 
these  metals.  As  a  result,  the  standard  Rice-Thomson  criterion  for  cleavage-like  crack 
propagation  [3]  pb/ys  >  7.5  —  10  is  fulfilled  for  these  metals  due  to  relatively  large  elastic 
moduli  in  comparison  with  the  decohesion  energy. 


In  order  to  identify  fundamental  characteristics  driving 


Table  1:  Cauchy  pressure  (C|2-CM)/C44,  surface  energy  for 
{111}  plane  and  Rice-Thomson  ratio  /y„b b/ys  for  select  refractory 
Ll2  alloys  in  comparison  with  NijAl  and  ductile  fee  metals. 


Metal  / 
alloy 

^12  ^44 

C44 

yUnb b  (Gpa) 

Ys  0/m2) 

iUnb^Ys 

Ni 

0.15 

68.5 

2.3 

7.6 

A1 

1.15 

24.8 

1.2 

5.9 

Ni3Al 

0.18 

57.5 

1.9 

7.8 

Ir 

-0.05 

198 

2.5 

24.5 

Rh 

-0.05 

133 

2.5 

14.6 

Ir3Nb 

-0.10 

173 

2.6 

18.4 

Ir3Zr 

0.05 

134 

2.4 

15.0 

Rh3Nb 

-0.03 

124 

2.3 

14.8 

Rh3Zr 

0.11 

87.6 

2.0 

12.2 

the  fracture  behavior  of  Ir-  and  Rh- 
based  alloys,  we  performed  ab  initio 
calculations  of  decohesion  energies 
and  elastic  moduli  in  Ir3Nb,  Ir3Zr, 
Rh3Nb,  and  Rh3Zr  alloys  by  means 
of  the  all-electron  full-potential 
linearized  augmented  plane  wave 
(FLAPW)  method  [4], 

The  surface  energy  ys  was 
determined  from  the  asymptotic 
value  of  the  ideal  cleavage 
decohesion  energy  =2yt, 

which  required  cleaving  an  infinite 
bulk  crystal  into  two  semi-infinite 
parts  and  the  creation  of  two  free 
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surfaces.  One  can  see  from  Table  1,  Ir-  and  Rh-based  alloys  (except  for  Rh3Zr)  as  well  as 
pure  Ir  and  Rh,  have  negative  Cauchy  pressures,  (C|2-C44)/C44.  This  feature  is  in  contrast 
with  the  Ni-based  alloys  and  confirms  the  pseudo-covalent  chemical  bonding  in  Ir3X  and 
Rh3X  intermetallics.  We  found  that  in  all  the  alloys  considered,  the  minimal  surface  energy 
corresponds  to  the  {111}  plane.  Though  the  shear  modulus  pnb  (where  n  is  the  normal  to  the 
{111}  slip  plane,  b  is  shear  direction  (110))  in  refractory  alloys  was  found  to  be  3  times 
larger  than  in  Ni3AI,  the  surface  energy,  ys,  is  only  approximately  25%  larger.  As  a  result, 
the  ratio  p„\,blys  for  Ir-  and  Rh-based  alloys  exceeds  significantly  the  Rice-Thompson  critical 
value.  Thus,  based  on  the  results  of  these  calculations,  we  predict  cleavage-like  crack 
propagation  in  these  alloys. 

2.  Superdislocation  structure  in  Rh-based  LI 2  alloys. 

Two-phase  y/y’  nickel-based  superalloys  remain  important  structural  materials  for  high- 
temperature  applications.  However,  their  potential  for  a  further  increase  of  the  temperature 
capability  of  these  materials  has  been  exhausted.  Recently,  a  new  class  of  superalloys  based 
on  platinum  group  metals  (PGM)  with  higher  melting  temperatures  and  superior 
environmental  properties  was  proposed  [5].  Among  these,  Pt-,  Ir-  and  Rh-based  alloys  attract 
the  most  attention.  Although  Rh-based  alloys  have  lower  melting  points  compared  to  Ir 
alloys,  they  possess  a  large  advantage  in  density.  Moreover,  it  has  been  known  that  the 
oxidation  resistance  of  Rh  is  superior  to  Ir  [6].  Although  there  are  few  doubts  as  to  the 
potential  of  the  PGA  for  ultra-high  temperature  applications,  further  progress  in  their 
development  is  hindered  due  to  an  insufficient  understanding  of  the  fundamental  factors  that 
control  their  mechanical  properties.  In  particular,  there  is  no  reliable  information  about 
dislocation  structure  and  mobility  in  Ir  and  Rh-based  alloys,  except  for  a  single  observation 
of  superdislocations  in  Ir3Nb  [7], 

Recently,  we  used  our  theoretical  combined  ab  initio  and  Peierls-Nabarro  (PN)  model 
approach  to  investigate  the  dislocation  properties  in  Ir3X  alloys  [8]  and  demonstrated  that 
features  of  the  dislocations  in  Ir-based  alloys  are  closely  connected  with  instability  in  relation 
to  Ll2— »  DO19  transformation.  Now  we  focus  on  Rh  group  alloys  and  apply  highly  accurate 
first-principles  calculations  of  the  generalized  stacking  fault  (GSF)  energetics  to  relate  the 
microscopic  characteristic  of  interatomic  interactions  with  properties  of  dislocation  in  Rh3X 
intermetallics  (where  X  =  Ti,  Zr,  Hf,  V,  Nb,  Ta). 


uA  uA 


Figure  1 :  Section  of  the  GSF  surface  in  (a)  lr3X  and  (b)  Rh3X  alloys  along 
(I  I2>  direction  corresponding  to  the  path  from  ideal  APB  to  SISF  positions. 


In  order  to  determine  the 
GSF  energies,  <I>(u)  which 
are  associated  with  a  rigid 
shift  of  one  part  of  the 
crystal  with  respect  to 
another  on  a  fault  vector  u 
in  the  slip  plane,  we 
carried  out  total  energy 
first-principles  calculations 
using  the  all-electron 
FLAPW  method  [4]  and  a 
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supercell  geometry  with  6  atomic  layers  and  homogeneous  periodic  boundary  conditions. 
For  RI13X  we  found  a  similar  tendency  as  was  observed  earlier  for  Ir3X  intermetallics  [8]  (cf.. 
Fig.  1).  All  alloys  can  be  divided  into  two  groups:  Rh  alloys  with  IVA  elements  (Hf,  Zr,  Ti) 
which  have  a  “normal”  ratio  of  the  superlattice  intrinsic  stacking  fault  (SISF),  y,  to  antiphase 
boundary  (APB)  energy  and  alloys  with  VA  elements  (Ta,  Nb,  V)  which  have  a  small 
y/^in  ratio.  The  decrease  of  the  SISF  energy  by  replacing  of  IVA  alloying  elements  with  VA 
has  an  electronic  origin  [8]  and  is  a  common  feature  of  both  Ir-  and  Rh-based  alloys. 

For  analyses  of  the  dislocation  properties,  we  employ  our  combined  approach  based  on  ab 
initio  y-surface  energetics  and  larger  scale  modeling  within  the  PN  model.  The  dislocation 
structure  is  shown  in  Fig.  2  as  a  “splitting  path”  together  with  the  y-surface  for  Rh3Ta,  Rh3Zr 
and  Rh3Nb.  We  found  that  superdislocations  with  type  II  core  structure  (Rear  splitting 
scheme  with  the  SISF  ribbon  between  partials)  are  strongly  preferred  energetically  in  Rh3Nb 
(Fig.  2(c)),  and  type  I  superdislocations  (with  APB  ribbon  between  partials)  are  predicted  for 
Rh3Ta  (Fig.  2(a))  and  Rh3Ti. 

The  Rh3X  intermetallics  are  distinguished  from  Ir3X  by  a  relatively  small  value  of  the  energy 
barrier  between  APB  and  SISF  (Fig.  1);  the  barrier  practically  disappears  for  Rh3Hf  and 
Rh3Zr.  As  a  result,  an  unusual  superdislocation  core  structure  was  found  in  Rh3Hf  and  Rh3Zr 
(Fig.  2(b)).  Because  of  the  very  shallow  APB  minima,  type  I  dislocations  become  unstable, 
and  SISF  bounded  type  II'  (with  different  sequence  of  Shockley  partials  [9])  is  preferable. 

Thus,  among  the  Rh3X  intermetallics  we  predict  that  type  I  dislocations  will  be  present  only 
in  Rh3Ti  and  Rh3Ta,  and  type  II  superdislocations  in  other  cases.  There  are  two  reasons  for 
the  appearance  of  SISF  bounded  superdislocations:  (i)  a  low  y/£m  ratio,  and  (ii)  an  unstable 
APB  energy;  these  conditions  result  in  the  appearance  of  type  II  and  type  IF 
superdislocations,  respectively. 


Figure  2:  Contour  plots  of  the  generalized  stacking  fault  energies  for  Ll2  (a)  Rh3Ta,  (b)  RhjZr  and  (c)  Rh3Nb; 
the  structures  of  the  superdislocations  are  shown  as  “splitting  paths”  (bold  curves). 

3.  Temperature  dependence  of  the  lattice  misfit  in  y/y'  superalloys.  Role  of  thermal 
expansion  and  composition  changes. 

The  magnitude  of  the  lattice  misfit  between  the  y  and  y'  phases,  S  =  2(ay  -  ay)/(  ay  +  ar) 
(where  ay,  ay  are  the  lattice  constants  for  the  corresponding  phases),  is  one  of  the  key 
parameters  which  determines  the  microstructure  formation,  stability  and  mechanical  behavior 
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of  modern  high-temperature  superalloys.  To  attain  a  stable  alloy  microstructure  ensuring  the 
high  creep  and  thermal  resistance,  the  misfit  parameter  should  have  an  optimal  value  at  the 
operational  temperature. 

Experimentally,  the  dependence  6(T)  in  the  Ni-based  alloys  exhibits  a  slight  decrease  in  the 
temperature  range  up  to  900°C  and  a  rapid  fall  for  T>  900°C  [10].  In  Ir-based  alloys,  the 
S(T)  shows  a  very  small  variation  in  the  temperature  region  up  to  1 500°C  [11]. 

The  relative  importance  of  two  contributions  to  the  temperature  dependence  of  6(T)  are  under 
intense  investigation,  namely:  (i)  the  difference  in  thermal  expansion  of  the  two  phases  [12], 
and  (ii)  the  redistribution  of  alloying  component  between  y  and  y'  with  the  increase  of 
temperature  [13]. 

We  explore  the  role  both  thermal  expansion  and  composition  changes  in  d(T)  for  the  Ni-AI 
and  Ir-Nb  y/y'  two-phase  alloys  based  on  ab  initio  total  energy  and  phonon  spectra 
calculations.  The  thermal  expansion  and  lattice  misfit  were  calculated  in  the  framework  of 
the  quasi-harmonic  approximation  on  the  basis  of  the  linear  response  theory  for  phonon 
spectra  in  Ni,  Ir,  Ni3AI,  Ir3Nb  using  the  pseudopotential  method  [14].  High  precision 
calculations  of  the  lattice  parameter  changes  due  to  AI(Nb)  substitutions  in  the  Ni(Ir)  matrix 
and  antisite  defects  in  the  y'  phase  were  performed  on  the  basis  of  total  energy  calculations 
using  the  highly  accurate  FLAPW  method  [4], 

The  results  of  calculations  of  the  thermal  expansion  coefficients,  a(T),  and  the  temperature 
dependence  of  the  bulk  moduli,  B(T),  for  Ni,  Ni3Al,  Ir,  I3Nb  (Fig.  3)  are  in  good  agreement 
with  available  experimental  data.  A  small  deviation  from  experiment  for  Ni  and  Ni3AI 
indicates  that  the  quasi-harmonic  approximation  becomes  less  accurate  at  high  temperatures. 
Our  results  also  clearly  demonstrate  the  important  contribution  of  the  ferromagnetism  of  Ni 
to  its  B(T)  behavior. 


(b,c,e,f)  for  Ni/NijAl  (a-c),  Ir/Ir3Nb  (d-f)-  The  symbols  mark  the  experimental  data. 
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Since  the  values  of  the  thermal  expansion  coefficients  a(T)  are  very  similar  between  the  y 
and  corresponding  y'  phases  (Fig.  3(a), (d)),  the  thermal  expansion  results  in  a  very  slow 
variation  of  the  lattice  misfit  d(T)  with  temperature  (Fig.  4(a), (b)).  Such  behavior  is  in 
agreement  with  experiment  for  lr/lr3Nb  (Fig.  4(b)),  but  drastically  different  from  the 
observed  S(T)  dependence  for  the  Ni/Ni3AI  alloys  (Fig.  4(a)).  The  latter  indicates  that  the 
redistribution  of  alloying  elements  between  phases  could  be  the  main  contributor  to  the 
lattice  misfit  dependence  at  high  temperatures.  The  concentration  dependence  of  the  lattice 
parameters  can  be  presented  in  the  form: 

ar{c)=a(*\\+cYpY),  ^(c)  =  40)( \+(C/  -0.25)4), 

where  cy,  cy  are  concentrations  of  the  redistributed  component,  and  /?y,  fiy-  are  the 
concentration  expansion  coefficients,  da/dc ,  in  the  y  and  y'  phases,  respectively.  The  values 
of/?  were  obtained  from  FLAPW  calculations  of  the  volume  changes  of  32-atom  supercells 
with  3.125%  of  the  redistributed  component,  and  are  as  follows:  fiy= 0.059,  fiy-=0.049  for 
Ni/Ni3AI  and  4=0.062,  4=0.04  for  Ir/Ir3Nb.  The  temperature  dependence  of  the 
concentrations  cy(T),  cy(T)  were  assumed  to  follow  the  solubility  lines  in  the  y  and  y'  phases 
on  the  phase  diagram  (curves  I,  2  in  Fig.  4(c), (d));  this  is  true  for  two-phase  systems  in 
thermal  equilibrium.  The  calculated  effect  of  the  composition  changes  is  presented  in  Fig. 
4(a).  It  is  seen  that  its  slope  is  very  close  to  the  experimental  values  at  high  temperatures.  At 
low  temperatures  (<  600  K)  where  the  thermal  diffusion  is  suppressed,  the  lattice  misfit 
follows  the  behavior  predicted  from  the  thermal  expansion  alone,  i.e.  it  increases  slightly. 
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Figure  4:  Contributions  of  the  thermal  expansion  and  composition  changes 
to  the  temperature  dependence  of  the  constrained  lattice  misfit  6C(T)= 
6(T)-3B/(3B+4p)  for  (a)  Ni/NijAl  and  (b)  Ir/lrjNb  (p  is  the  shear  modulus). 
Rectangles  correspond  to  the  experimental  data;  (c)  Ni-rich  and  (d)  Ir-rich 
parts  of  equilibrium  phase  diagrams  for  Ni-AI  and  Ir-Nb  alloys,  respectively. 


between  the  alloy  phase  diagram  and  the  temperature  dependence  d(T). 


Thus,  the  redistribution  of 
the  major  alloy  components 
(Al  into  Ni  and  Ni  into 
Ni3AI)  gives  the  main 

contribution  to  the  tempe¬ 
rature  dependence  of  the 
lattice  misfit  and  the 
behavior  of  d(T)  for 

Ni/Ni3AI  at  high  tempe¬ 
ratures.  For  the  Ir/Ir3Nb 
system,  the  alloy  compo¬ 
nent  redistribution  will  start 
to  contribute  to  the  lattice 
misfit  only  at  extremely 
high  temperatures  (>2000 
K).  The  amplitude  of  these 
contributions  can  be 

determined  by  considering 
the  shape  of  the  y  -  y'  gap  on 
phase  diagrams.  The  last 
conclusion  is  important  for 
alloy  design  as  it  allows  one 
to  establish  a  simple  relation 
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4.  The  electronic  origin  of  solid  solution  softening  in  bcc  molybdenum  alloys. 

The  addition  of  elements  which  improve  the  room  temperature  strength  and  high  temperature 
creep  resistance  of  materials,  is  important  in  alloy  design.  Besides  the  solid  solution 
hardening  effect  (SSH),  which  is  connected  with  the  pinning  action  of  impurities  on 
dislocations,  in  some  cases  the  alloying  can  lead  to  the  opposite  phenomenon,  the  decrease  of 
hardness,  which  is  known  as  solid  solution  softening  (SSS).  The  SSS  effect  is  well  studied  in 
bcc  transition  metal  alloys  and  attracts  remarkable  attention  as  a  possible  way  to  improve  the 
low  temperature  ductility  of  the  refractory  group  VIA  bcc  metals  (Cr,  Mo,  W)  [15]. 

It  is  commonly  accepted,  that  the  intrinsic  solution  softening  is  connected  with  direct 
impurity-dislocation  interactions  and  results  from  the  enhanced  kink  nucleation  at  point 
defects  [16].  However,  the  microscopic  mechanism  of  this  phenomenon  still  remains  a  matter 
of  discussion.  In  particular,  the  changes  in  hardening  rate,  dH/dc,  for  Mo  alloys  correlate 
with  atomic  size  misfit  only  for  W,  Ta,  Hf  leading  to  SSH  and  not  for  Pt,  Ir,  Os,  Re,  which 
promote  softening  for  low  solute  concentration  and  low  temperatures  [  1 5],[  1 7].  Moreover, 
for  binary  and  ternary  Mo  alloys  there  is  a  strong  correlation  between  softening  rates  and  the 
number  of  extra  electrons  brought  by  solute  elements  [17].  These  results  clearly  indicate  the 
importance  of  the  “chemical”  mechanism  of  impurity-dislocation  interactions  in  the  SSS. 


Here  we  present  an  investigation  of  the  microscopic  mechanisms  of  SSS  in  binary  Mo  alloys 
on  the  basis  of  ab  initio  electronic  structure  calculations.  Our  calculations  are  inspired  by  the 
observation  that  SSS  results  from  the  addition  of  certain  alloying  elements  which  decrease 
the  stacking  fault  energy  (SFE)  [1 8]  -  one  of  the  key  parameters  which  control  the  dislocation 
structure.  We  calculated  the  GSF  energies  in  molybdenum  alloying  with  5r/and  3d  transition 
metal  elements  and  showed  that  the  SSS  in  Mo  alloys  can  be  explained  by  the  decrease  of  the 
energy  barrier  for  kink  nucleation  due  to  the  local  lowering  the  SFE  at  the  substitutional 
atoms.  The  total  energy  calculations  were  carried  out  by  the  first-principles  full-potential 
linear  muffin-tin  orbital  (FLMTO)  method  [19]. 


The  GSF  energies 
estimated  for  (100)  planar 
shifts  along  (111)  (parallel 
to  the  Burgers  vector)  are 
shown  in  Fig.  5  for  bcc  Mo 
and  Mo  alloys.  All  5 d  TM 
additions  with  an  extra 
electron  concentration  lead 
to  the  reduction  of  the  GSF 
energy  in  comparison  with 
unalloyed  molybdenum. 
The  energy  of  the  unstable 
stacking  fault,  yus  (cor¬ 
responding  to  the  maximum  of  the  GSF  energy  in  Fig.  5),  decreases  monotonically  with  the 
atomic  number  of  the  solute  atoms  (cf.  Table  2).  The  additions  that  reduce  the  electron 
concentration  (Ta,  Hf)  as  well  as  an  isoelectonic  impurity  (W)  lead  to  an  increase  in  the  GSF 


Figure  5:  GSF  energies  for  <11 1X1 10)  shear  in  Mo  and  Mo  alloys  with 
(a)  5 d  and  (b)  3 d  elements.  The  x  axes  are  given  in  units  of  the  Burgers 
vector  1  /2<  111). 


7 


Table  2:  The  ratio  of  stacking  fault  energies  y, Jy°„  (/,„  corresponds  to  unalloyed  Mo),  the  ratio  of  the  lj/eg 
occupations  for  addition  atoms,  the  ratio  of  ionic  radii  R,/RMo  and  initial  hardening  rate  [5],  dH/dc ,  for  binary 
molybdenum  alloys. 


Hf 

Ta 

W 

Re 

Os 

Ir 

Pt 

y  us^i  us 

1.49 

1.09 

1.02 

0.86 

0.63 

0.46 

0.38 

t2/eg 

1.86 

1.85 

1.86 

1.81 

1.73 

1.58 

1.41 

R/Rmo 

1.13 

1.05 

1.01 

0.99 

0.96 

0.97 

0.99 

dH/dc  \6l 

10 

5 

~0 

-5 

-35 

-70 

-120 

energy  and  the  additions  with  excess  electrons  (Re,  Os,  Ir  and  Pt)  lower  its  value.  The 
decrease  in  GSF  energy  obtained  in  the  row  of  solutes  from  Re  to  Pt  strongly  correlates  with 
the  experimentally  observed  lowering  of  the  initial  hardening  rate  [17]  (i.e.  dH/dc  for  C 
approaching  zero)  at  77  K  (Table  I).  Thus,  the  calculated  results  point  out  that  the  SSS 
mechanism  in  Mo  alloys  should  be  connected  with  the  decrease  in  GSF  energy,  similar  to 
that  established  in  [18]  for ycc-based  alloys.  Similar  trends  were  obtained  for  3r/TM  solutes, 
but  in  this  case  all  elements  including  Ti  and  V  lead  to  a  decrease  of  yu.  Hence,  based  on 
these  results  one  can  also  predict  the  SSS  effect  for  3 d  solutes  with  extra  electrons. 

To  make  clear  the  relation  between  the  chemical  bonding  and  the  GSF  energy  characterizing 
the  lattice  resistance  to  (1  1  1)(1 10)  shear  deformation,  we  analyzed  the  changes  in  the  orbital 
occupations  for  additions  in  Mo  binary  alloys.  For  the  bcc  structure,  the  l2g  and  eg  orbitals  are 
responsible  for  the  directional  bonds  between  the  nearest  neighbor  atoms  ( nn ,  along  the  (111) 
direction)  and  second  (next  nearest)  neighbors  ( nnn ,  along  (100)),  respectively.  We  found 
that  the  local  changes  in  the  occupation  anisotropy  correlate  with  the  GSF  energy  (cf.,  Table 
1)  and  point  out  that  the  decrease  in  the  shear  resistance  is  connected  with  the  weakening  of 
the  nn  bonds  relative  to  nnn  near  the  solute  atom. 


Thus,  the  calculated  trends  in  the  GSF  energies  confirm  the  suppositions  based  on  the 
experimental  investigations  [17]  concerning  the  electronic  nature  of  the  SSS.  This  conclusion 
is  not  trivial  because  there  is  no  simple  relation  between  the  GSF  energy  and  kink  formation 
in  bcc  metals.  As  established  in  recent  calculations  (see  [20]  and  references  therein),  the 
screw  dislocation  in  Mo  has  a  nondegenerate  core,  spread  symmetrically  into  the  three  (110) 
planes.  We  suggest  that  a  solute  atom  locally  disturbs  the  three-fold  symmetry  of  the 
dislocation  core  due  to  the  lowering  of  the  GSF  energy  for  one  of  the  three  (110)  planes.  The 
barrier  for  double  kink  nucleation  can  be  roughly  estimated  from  the  energy  difference  A Ej  3 
between  a  plane  core  with  a  solute  atom  and  the  three-fold  core  in  the  matrix  [21], 


A£2  3=  —  In 
2-3  4/z- 


i-2y, 

v° 

•  us 


\  f  US  J 

Both  3d  and  5 d  transition  metals  lead  to  a  monotonous  decrease  in  A Ei  3  and  the  energy 
barrier  for  double  kink  nucleation  with  atomic  number  for  additions  with  extra  electrons. 
Thus,  we  demonstrated  that  the  tendencies  in  SSS  can  be  understood  in  terms  of  the  effect  of 
solute  atoms  on  the  GSF  energy. 
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5.  Effect  of  additions  on  dislocation  structure  and  mobility  in  bcc  Mo.  Atomic  row  model 
with  ab  initio  parameterization. 

A  strong  correlation  between  the  softening  behavior,  unstable  stacking  fault  energy  and 
number  of  valence  electrons  introduced  by  alloying  additions  in  bcc  Mo  alloys  discussed 
above  clearly  points  out  an  electronic  origin  of  the  solid  solution  softening  and  emphasizes 
the  importance  of  the  changes  in  dislocation  structure  due  to  solutes. 


In  bcc  metals,  plastic  deformation  is  realized  by  the  motion  of  screw  dislocations  with  a  star- 
like  nonplanar  core,  and  being  a  planar  defect  the  stacking  fault  does  not  characterize  the 
dislocation  structure  quite  accurately.  To  investigate  the  effect  of  additions  on  the  core 
structure,  we  employ  a  combined  approach  that  includes  the  atomic  row  model  for  a  screw 
dislocation  in  bcc  metals  [22]  with  ab  initio  parametrization  of  the  interatomic  row 
potentials,  O(w).  The  potential  <D(m)  is  assumed  to  be  a  function  only  of  the  relative  atomic 
positions,  u,  and  is  composed  of  the  first  and  second  terms  in  the  Fourier  expansion: 


.  .  .  Oo  (  2 nu  \nu  1  „ 

q>(w)  = - 51 —  cos - or  cos - +  -(  -or) 

(l-4«)v  b  b  2K  ' 


(1) 


where  b  is  the  Burgers  vector.  The  energy  variation  under  the  atomic  row  shift  along  (111), 

A  Em  («)  =  3(4>(«  - 1/3)  +  0(1/  + 1/3)),  (2) 

was  calculated  using  the  first-principles  full-potential  linear  muffin-tin  orbital  (FLMTO) 
method  [19]  and  the  parameters  <D0  and  or  found  from  Eqs.  (1)  and  (2).  The  dislocation  core 
structure  is  determined  by  the  relaxation  of  each  atomic  row  position,  w„  with  appropriate 
boundary  conditions  to  its  equilibrium  value  to  provide  a  minimum  of  the  total  energy. 


0  0.2  0.4  0.6  0.8  1 

displacement  u/t 


The  calculated  atomic  row  displacement  (ARD) 
energy,  A Enne,  (Fig.  6)  decreases  for  the  solutes  with 
extra  electrons  (Re,  Os,  Ir,  Pt)  and  increases  for 
additions  that  reduce  the  electron  concentration  (Hf, 
Ta  and  W).  The  changes  in  atomic  row  shift 
resistance  correlate  with  the  variation  of  yus  [23]  and 
with  the  observable  trend  in  SSS  [24],  The  effect  of 
additions  is  not  localized  in  the  atomic  row  with  the 
solute,  but  covers  the  region  of  the  nearest  atoms 
around  the  solutes. 


Figure  6:  Atomic  row  displacement  energies.  We  demonstrated  that  the  atomic  row  model  with 
A£„w,  for  a  shift  of  the  atomic  row  with  ab  initio  parametrization  of  the  interatomic  row 
solute  along  (111)  versus  displacement  in  potential  gives  the  correct  description  of  the 
units  of  the  Burgers  vector  1/2(111).  dislocation  structure.  In  particular,  we  found  a 

compact  isotropic  core  structure  in  Mo  that  is  in 
agreement  with  the  results  of  direct  first  principles  calculations  [25].  We  modeled  the  effect 
of  the  solute  by  the  relevant  change  in  the  interatomic  row  potential  and  found  that  the  solute 
in  the  “easy”  dislocation  core  disturbs  the  three-fold  symmetry  due  to  (i)  local  distortions  in 
Re  and  Os  or  (ii)  a  transformation  to  the  split  (Fig.  7a, b)  or  spread  (Fig.  7c)  core  structure  for 
Ir  and  Pt  additions.  The  softener  additions  (with  extra  electrons)  attract  the  dislocation  in 
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Figure  7:  Effect  of  Pt  additions  (shown  by  filled  circle)  on  the  “easy”  dislocation  core  structure:  (a)  slightly 
distorted  core  due  to  Pt  in  position  3  far  from  the  dislocation  core,  (b)  split  core  structure  due  to  solute  in 
position  2,  (c)  spread  core  structure  with  solute  in  the  dislocation  core  (position  I ). 

contrast  with  hardeners  which  have  repulsive  solute- 
dislocation  interactions.  As  a  result,  softener  additions 
will  promote  double  kink  nucleation  because  they 
produce  an  extra  driving  force  (Fig.  8)  for  bowing  the 
dislocation  line. 

Thus,  the  atomic  row  model  predicts,  in  agreement  with 
experiment,  the  enhancement  of  double  kink  nucleation 
due  to  softener  solutes.  Because  ARD  energies 
gradually  change  with  the  band  filling  in  alloys  for 
softener  additions,  our  Findings  allow  one  to  explain  the 
observed  correlation  between  SSS  and  the  atomic 
number  of  the  solute  [24], 

6.  Solid  solution  softening  and  hardening  in  the  Group  V  and  Group  VI  refractory  hcc 
transition  metal  alloys.  Atomic  row  model  with  ab  initio  parameterization. 

In  contrast  with  the  Group  VIA  alloys,  the  softening  in  the  Group  VA  (V,  Nb,  Ta)  bcc  metal 
alloys  was  traditionally  considered  as  an  extrinsic  phenomenon  originating  from  the 
scavenging  of  interstitial  impurities  by  the  solutes  [26].  Recently,  however,  the  Peierls 
mechanism  was  suggested  as  responsible  for  hardening  rates  in  Ta-W  alloys  at  low 
temperatures  [27],[28].  As  result,  the  question  about  the  microscopic  mechanism  for  the 
hardening/softening  behavior  in  the  Group  VA  bcc  refractory  metals  is  still  open  to 
discussion.  Here  we  present  the  results  of  comparative  investigations  of  the  solution 
softening  in  the  Groups  VA  (Nb,  Ta)  and  VIA  (Mo,  W)  metal  alloys.  Direct  investigations  of 
the  solute-dislocation  interaction  by  first  principles  electronic  structure  calculations  are  still  a 
difficult  problem.  To  make  clear  the  effect  of  additions  on  the  dislocation  core  structure  and 
mobility,  we  employ  our  combined  approach  [23],[29],  which  includes  the  atomic  row  model 
(ARM)  for  a  screw  dislocation  in  bcc  metals  with  the  atomic  row  displacement  energy, 
A Erow,  calculated  using  the  ab  initio  full-potential  linear  muffin-tin  orbital  (FLMTO)  method 

[19]. 


A  n 


Figure  8:  The  contribution  to  the  driving 
force  for  double  kink  nucleation,  A Ejf,  for 
Pt,  Re,  and  Hf  additions. 


The  dependence  of  A Erm(u)  on  the  relative  atomic  row  displacement  u  is  shown  in  Fig.  9. 
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The  value  of  A Erow  decreases  gradually  with  the  number  of  extra  valence  electrons  of  solutes 
(Re,  Os,  Ir,  Pt)  in  both  Mo  and  W  (Fig.  9(a,b))  and  these  trends  correlate  with  the 
experimental  initial  hardening  rates.  For  the  Group  VA  metals,  we  obtained  the  opposite 
result:  the  Ad,  5 d  additions  with  extra  valence  electrons  increase  the  atomic  row  shift 
resistance  and  the  additions  with  a  lack  of  electrons  reduce  A Erow(u).  This  difference 
originates  from  the  peculiarities  of  the  electronic  structure  and  chemical  bonding,  when  the 
Fermi  level  is  located  in  the  region  of  the  bonding  states  (Group  VA)  or  in  the  pseudogap 
between  bonding  and  antibonding  states  (Group  VIA). 
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Using  the  atomic  row 
model  with  ab  initio 
parameterization,  we  found 
that  the  screw  dislocations 
in  pure  Mo,  W  as  well  as  in 
Nb  and  Ta  metals  have 
isotropic  compact  cores. 
The  driving  forces  for  the 
double  kink  nucleation  - 
determined  as  total  energy 
changes  when  a  dislocation 
moves  from  a  far  position 
to  the  solute  -  are  presented 
in  Fig.  10.  The  softener 
additions  which  decrease 
the  ARD  energy  have  an 
attractive  interaction  with 
the  dislocation,  in  contrast 
with  alloying  elements  that 
increase  A Erow(u)  and  have 
a  repulsive  interaction  with 

dislocations.  For  the  Group  VA  alloys,  the  effect  of  solutes  on  the  driving  forces  is  weaker 
and  does  not  demonstrate  the  monotonic  dependence  on  the  number  of  valence  electrons 
(Fig.  10).  The  additions  with  extra  electrons  in  the  Group  VIA  metals  enhance  the  double 
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Figure  9:  ARD  energies,  A Em„,  for  a  shift  of  atomic  row  with  solute  along 
(100)  in  (a)  Mo,  (b)  W,  (c)  Nb  and  (d)  Ta  alloys. 


kink  nucleation  and  result  in  softening,  but  they  lead  to 


0  02 


001 


e 


-0.01 


-0.02 


-2 


A  n 


2  4 

A  n 

Figure  10:  Driving  force  for  the  kink  nucleation  as  a  function  of  the 
difference  in  the  number  of  valence  electrons  for  (a)  the  Group  VIA  alloys 
( 1  -  Mo,  2  -  W)  and  (b)  Group  VA  alloys  (I  -  Nb,  2  -  Mo). 


a  repulsive  interaction  with 
dislocations  and  cannot  be 
intrinsic  softeners  in  the  Nb 
and  Ta  alloys.  By  contrast, 
the  solutes  that  lack 
electrons  suppress  the  kink 
formation  in  the  Group  VIA 
metals,  and  promote  it  in 
the  Group  VA  metals.  In 
agreement  with  experiment 

[30],  a  strong  strengthening 
effect  in  Nb  and  Ta  was 
found  for  the  Re  additions. 
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in  contrast  with  the  Group  VIA  alloys  where  Re  results  in  softening  [15]. 

Thus,  our  results  allow  one  to  conclude  that  the  effect  of  alloying  with  4 d  and  5 d  transition 
elements  on  dislocation  structure  and  mobility  is  fundamentally  different  for  the  bcc 
refractory  metals  of  Groups  VIA  and  VA.  The  additions  with  extra  valence  electrons,  which 
give  rise  to  SSS  in  the  Group  VIA  metals,  cannot  lead  to  intrinsic  softening  in  the  Group  VA 
metals.  Oppositely,  the  additions  with  a  lack  of  electrons  (leading  to  strong  hardening  in 
Group  VI  metals)  may  result  in  intrinsic  softening  in  the  Group  V  metals,  which  is  controlled 
by  the  double  kink  nucleation  mechanism.  We  believe  that  the  softening  observed  in  Group 
V  metals  alloyed  with  Mo,  W,  and  Re  is  governed  instead  by  the  extrinsic  mechanisms  -  for 
example,  by  scavenging  of  solute  interstitials.  For  the  first  time,  we  demonstrated  that  both 
softening  and  hardening  effects  in  the  bcc  refractory  metals  are  determined  by  the  electronic 
structure  rather  than  by  atomic  size  or  shear  modulus  misfits  which  are  traditionally 
considered  as  the  main  parameters  controlling  the  hardening/softening. 

7.  Electronic  mechanism  of  the  stabilization  of  A15  phase  in  Cr  and  Mo  alloys  by  light 
impurities. 

Alloys  based  on  the  group  VI-A  refractory  bcc  metals  (Cr,  Mo,  W)  are  attractive  materials 
for  high-temperature  applications.  Among  them,  chromium  alloys  have  advantages  due  to 
low  density,  high  creep  and  oxidation  resistance  [31],  and  molybdenum  alloys  attract  special 
attention  due  to  their  potential  for  aerospace  applications.  Unfortunately,  the  lack  of  ductility 
at  low  temperatures  has  severely  limited  the  use  of  bcc  refractory  metals,  whose  room 
temperature  brittleness  is  caused  by  extrinsic  factors  -  extremely  low  solubility  of  interstitial 
impurities  (C,  N,  O)  resulting  in  the  formation  of  brittle  carbides,  nitrides  and  oxides. 

Our  recent  ab  initio  calculations  [32]  showed  that 
small  Cr-Re  clusters  with  TCP  A15  structure  may 
become  energetically  favored  by  the  trapping  of 
carbon  impurities.  This  finding  supports  a  proposed 
mechanism  [33]  of  the  “rhenium  ductilizing  effect” 
which  connects  the  increase  of  the  solubility  of 
impurities  with  the  appearance  of  TCP  particles. 
These  particles  serve  as  efficient  “attractors”  of 
interstitials  and  prevent  the  formation  of  brittle 
phases  and  grain  boundary  segregation. 

Further,  we  evaluated  the  role  of  light  impurities  (O, 
N,  C)  on  the  stabilization  of  the  A15  phases  of  Cr 
and  Mo  on  the  basis  of  first-principles  FLAPW  [4] 
calculations.  The  A15  phases  of  pure  Cr  and  Mo  are 
metastable,  higher  in  energy  than  the  bcc  phase  by 
only  70  and  84  meV/atom  for  Cr  and  Mo, 
respectively.  For  the  A15  Cr  and  Mo  phases  with 
oxygen  impurities,  both  substitutional  and  interstitial 
positions  were  investigated.  In  all  cases,  full  atomic 
relaxations  were  performed,  and  lowest  energy 
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Figure  11:  Electronic  densities  of  states  for 
A15  Cr  with  interstitial  O. 
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configurations  were  obtained.  Substitutional  structures  were  found  to  be  very  unstable,  while 
the  formation  of  the  interstitial  A 15-0  phase  is  energetically  favored.  As  seen  from  Fig.  1 1, 
the  stabilizing  effect  of  interstitial  oxygen  in  A15  Cr  is  a  result  of  strong  hybridization  of 
02 p  and  Cr3 d  states  over  the  whole  bonding  energy  range,  which  also  manifests  itself  by  the 
opening  of  a  sharp  and  deep  DOS  minimum  exactly  at  the  Fermi  energy.  For  the  A15  Mo 
phase,  a  similarly  strong  hybridization  of  02p  and  Mo4J  states  takes  place. 


Table  3:  Lattice  parameters  of  AI5  Cr  and  Mo  phases  with 
interstitial  impurities  and  their  stability  with  respect  to  bcc 
Cr(Mo)  with  interstitial  impurities. 


The  formation  energies  were 
calculated  for  A  15-type  structures  with 
O,  N  and  C  interstitials.  Table  3 
presents  calculated  lattice  parameters 
and  tetragonal ity  of  the  structures; 
these  compare  favorably  with  available 
experiments  on  metastable  A15  phases 
observed  in  thin  films  [34],  In  bulk  Cr, 
these  structures  will  compete  with 
interstitials  in  the  bcc  phase.  The 
calculated  energy  difference,  A E  = 
£(A15)  -  £(bcc)  (Table  3),  shows  that 
the  configurations  with  interstitials  in 
the  bcc  phase  are  still  more  energetically  favored  than  in  the  A15  phase,  and  oxygen 
interstitials  are  the  most  preferred  for  the  formation  of  the  A15  phase.  Thus,  A15  Cr  and  Mo 
are  impurity-locked  metastable  phases. 


Structure 

a  (A) 

c/a 

£(A15)-£|bcc) 

(eV/atom) 

A15Cr+Ojn, 

4.593 

1.058 

0.05 

A15Cr+Nj„, 

4.571 

1.063 

0.14 

A 1 5Cr+Cin, 

4.551 

1.080 

0.13 

A15Mo+Oin, 

5.069 

1.038 

0.08 

A15Mo+Ninl 

5.058 

1.039 

0.17 

A15Mo+Cinl 

5.054 

1.045 

0.20 

The  effect  of  alloying  Cr  and  Mo  as  a  means  to 
facilitate  the  formation  of  the  A15  phase  was  then 
investigated.  The  following  alloying  elements  were 
chosen  for  Cr  alloying:  Re,  a  known  ductilizer,  and  Fe 
and  Ni,  which  were  conceived  as  possible  Re 
alternatives.  The  calculated  formation  energies  for 
Cr3M  (M  =  Re,  Fe,  Ni)  compounds  in  the  A15 
structure,  shown  in  Fig.  12  as  a  dependence  of  the 
formation  energy  on  the  volume,  indicate  that,  while 
these  phases  remain  metastable,  they  have  lower 
energy  compared  to  the  A15  phase  of  Cr,  with  Cr3Re 


Figure  12:  Formation  energy  for  AI5  Cr3M 
phases  in  comparison  with  bcc  and  A15  Cr. 


showing  the  largest  energy  decrease.  For  Mo,  only 
the  effect  of  alloying  with  Re  was  investigated.  The 
calculations  of  Mo3Re  in  the  A15  structure  indicate 
that,  similarly  to  Cr,  it  has  lower  formation  energy  (+0.052  eV/atom)  compared  to  the  A15 
phase  of  Mo  (+0.084  eV/atom).  Thus,  the  alloying  with  Re  increases  the  stability  of  the  A 15 
phase  of  Mo. 


The  cumulative  effect  of  both  alloying  and  interstitial  oxygen  was  investigated  further.  In 
addition  to  the  formation  energies,  the  energy  gain  resulting  from  oxygen  interstitials  was 
also  calculated  as:  Egam  =  £’,£>,(Cr(Mo)3M)  +  Vi Elol  (O2)  -  Eto,  (Cr(Mo)3M+Ojnt).  As  seen  from 
Table  4,  the  energy  gain  is  highest  for  Cr3Re  and  Mo3Re,  where  it  is  2  and  6  times  higher 


13 


Table  4:  Energy  gain  due  to  O  interstitials  and  than  for  Cr3Ni  and  Cr3Fe;  respectively.  This  is  a 

result  of  both  strong  O-Re  bonding  and  lattice 
expansion  due  to  Re  alloying;  the  latter  creates  a 
larger  interstitial  pore  that  accommodates  the 
oxygen  atom  better  leading  to  a  smaller 
structural  deformation.  As  a  result,  the 
formation  energy  for  A15  Cr3Re  and  Mo3Re 
with  O  interstitials  becomes  largely  negative, 
confirming  the  stability  of  these  structures. 


Finally,  the  formation  energy  for  the 
Cr(Mo)3Re+Om,  phase  was  compared  with  that  of  the  O  interstitial  in  the  bcc  Cr(Mo)-25% 
Re  alloy,  and  the  preference  energies  Epref=  A15  Cr(Mo)3Re+Oint)  -  Elol( bcc  Cr(Mo)- 
25%Re+Omt)  were  calculated.  The  preference  energies  of  -0.064  eV/atom  and  -0.057 
eV/atom  were  obtained  for  Cr3Re+Omt and  Mo3Re+Ojm,  respectively.  Thus,  based  on  precise 
first-principles  total  energy  calculations,  we  demonstrate  that  Re  alloying  in  the  presence  of 
light  elements  may  result  in  the  formation  of  stable  Cr(Mo)3Re+interstitial  precipitate 
particles  with  A  15-like  structure.  These  particles  can  act  as  scavengers  of  light  impurities, 
increase  impurity  solubility  and  thus  directly  improve  the  ductility  of  Mo  alloys. 

8.  Fundamental  parameters  of  the  brittle-ductile  behavior  of  B2  intermetallics. 

The  transition  metal  (TM)  aluminides,  NiAl,  CoAl,  and  FeAl,  attract  considerable  attention 
as  potential  structural  materials  due  to  high  strength,  environmental  resistance  and  low 
density;  however,  they  show  brittleness  at  ambient  temperatures.  The  low  ductility  of  TM 
aluminides  comes  from  intrinsic  reasons  and  is  connected  with  their  susceptibility  to 
cleavage-like  crack  propagation  [35]  and  an  insufficient  number  of  slip  systems.  Recently,  a 
whole  new  class  of  ductile  B2-type  intermetallics  compounds,  based  on  rare-earth  elements 
(such  as  YCu,  YAg,  and  YZn)  has  been  discovered  [36].  In  addition,  high  ductility  was 
found  in  CoZr  [37],  which  belongs  to  a  unique  group  of  Co-based  B2  CoX  intermetallics  (X 
=  Ti,  Zr,  Hf)  exhibiting  a  positive  temperature  dependence  of  the  yield  stress  [38]. 

To  identify  fundamental 
characteristics  driving  the  fracture 
behavior  of  B2  intermetallics,  we 
calculated  the  decohesion 
energies  and  elastic  moduli  in  TM 
aluminides  and  Co-  and  Y-based 
alloys  (Table  5)  by  means  of  first- 
principles  total  energy  FLAPW 
[4]  calculations.  The  surface 
energy,  ys,  was  determined  from 
the  asymptotic  value  of  the  ideal 
cleavage  decohesion  energy 
ET  =  2 y„  required  to  cleave  an 
infinite  bulk  crystal  into  two 


Table  5:  Elastic  modulus  C44,  ratio  of  C44to  bulk  modulus  B , 
surface  energy  for  { 1 10}  plane  ys ,  and  Rice-Thomson  ratio  fib!yx 
for  TM  aluminide,  Co-  and  Y-based  B2  intermetallics. 


Alloy 

C44  (GPa) 

C44/B 

MJ/m2) 

pb/ys 

NiAl 

153 

0.85 

1.64 

27.6 

CoAl 

145 

0.81 

2.23 

18.8 

FeAl 

120 

0.76 

2.39 

14.4 

CoTi 

63 

0.35 

2.00 

9.4 

CoZr 

43 

0.28 

1.77 

7.6 

CoHf 

56 

0.35 

1.89 

9.3 

YCu 

33 

0.49 

1.15 

10.5 

YAg 

34 

0.52 

1.01 

12.2 

YZn 

40 

0.69 

1.05 

14.3 

formation  energy  of  A 15  Cr(Mo)3M+Oint. 


System 

Eform 

(eV/atom) 

Egain 

(eV/oxygen) 

Cr3Fe+Oin, 

+0.023 

0.24 

Cr3Ni+Oin, 

-0.026 

0.71 

Cr3Re+Oin, 

-0.127 

1.38 

Mo3Re+Oint 

-0.143 

1.56 
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semi-infinite  parts  and 
the  creation  of  two  free 
surfaces  (see  Fig.  13). 


It  is  seen  from  Table  5 
that  Co-  and  Y-based 
alloys  are  distinct  from 
TM  aluminides  by  the 
low  C44  shear  modulus. 
The  small  C44/B  ratio 
indicates  possible  lattice 
instability  with  respect 
to  martensitic  transformation.  Indeed,  YCu  and  CoTi  alloys  were  found  to  undergo  at  low 
temperatures  the  B2— >B27  and  B2— »B19  transformation,  respectively.  According  to  the 
Rice-Thomson  criterion  [3],  cleavage-like  (brittle)  crack  propagation  is  expected  if  the 
condition  fi b/ys  >  7.5-10  is  fulfilled.  (Here,  b  is  the  Burgers  vector,  fi  the  shear  modulus;  n  = 
C44  for  (100){  1 10}  shear).  Even  though  the  shear  moduli  of  Co-  and  Y-based  alloys  are  3-5 
times  lower  compared  to  TM  aluminides,  the  difference  in  the  surface  energy  ys  is  much  less. 
As  a  result,  the  ratio  nb/ys  is  significantly  higher  for  TM  aluminides  and  exceeds  the  Rice- 
Thompson’s  critical  value  resulting  in  cleavage-like  crack  propagation.  At  the  same  time,  we 
predict  ductile  behavior  for  Co-based  alloys  and  YCu  and  YAg,  due  to  the  plastic  relaxation 
processes  in  a  crack  tip,  and  low  fracture  toughness  for  YZn;  these  findings  agree  with 
experiment.  Thus,  the  observed  ductility  of  Co-  and  Y-based  B2  intermetallics  has  intrinsic 
origins  and  is  connected  with  the  low  shear  resistance  combined  with  the  relatively  high 
decohesion  energy. 

9.  Dislocation  structure  and  mobility  in  B2  CoZr:  combined  ab  initio-Peierls-Nabarro 
analysis. 


Figure  13:  Decohesion  energy,  Ect  (J/m2),  as  a  function  of  the  separation 
distance  u  (in  units  of  the  {110}  interplane  distance  dm)  in  (a)  TM  aluminides 
and  (b)  Co-  and  Y-  based  alloys. 


As  mentioned  above,  Co-based  intermetallics  satisfy  the  Rice-Thomson  criteria  for  ductile 
behavior,  and  remarkable  plasticity  was  observed  recently  in  CoZr  [37],  In  addition,  these 
alloys  are  unique  among  B2  intermetallics  due  to  a  positive  yield  stress  dependence  on 
temperature  [38],  despite  the  fact  that  only  (100)  dislocations  operate  in  plastic  deformation. 
To  understand  the  origin  of  the  unusual  mechanical  behavior  of  Co-based  alloys,  an 
investigation  of  the  dislocation  properties  is  necessary.  Experimentally,  there  is  only  one 
HREM  observation  of  dislocations  in  CoTi  [39], 


Table  6:  Relaxed  unstable  stacking  fault  energy 
Y"!"'  f°r  <100){01 1}  shear,  ratio  y,  /  y^lax  and 
width  of  dislocation  core  d,  for  CoZr  and  NiAI. 


Alloy 

r:,ax(  J/m2) 

yJyT 

d/a 

CoZr 

0.19 

9.8 

6.6 

NiAI 

1.36 

1.2 

2.4 

To  investigate  the  structure  and  mobility  of  (100) 
dislocations  in  B2  Co-based  intermetallics,  we 
employ  our  approach  which  combines  the 
modified  Peierls-Nabarro  model  with  ab  initio 
generalized  stacking  fault  (GSF)  energetics.  We 
found  that  GSF  energies  for  {001}  and  {110} 
planes  in  CoZr  are  significantly  smaller 
compared  with  NiAI,  and  the  lowest  shear 
resistance  corresponds  to  the  {110}  slip  plane 
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(previously,  we  reported  similar  behavior  for  CoTi).  Most  striking  is  a  dramatic  decrease  of 
the  GSF  energy  for  (100){0I  I }  shear  due  to  atomic  relaxation  in  CoZr,  in  contrast  with  NiAl 
where  the  GSF  energy  remains  high  after  relaxation  (Fig.  14).  In  particular,  the  unstable 
stacking  fault  energy,  yus,  decreases  from  1.07  to  0.18  J/nr  in  CoZr  for  the  {110}  plane  and 
from  1 .95  to  0.87  J/m2  for  the  {001 }  plane.  This  results  in  a  high  value  of  the  ys/yus  ratio  (see 
Table  6)  that  indicates  ductile  behavior  of  CoZr  according  to  a  criterion  proposed  in  [40], 

Taking  unto  account  the 
GSF  relaxation,  we 
predict  the  splitting  of 
the  <100){01 1 > 

dislocation  structure  in 
CoZr  to  occur  according 
to  the  following  scheme: 
b  =  b/2  +  b/2  (Fig.  1 5a). 
Additionally,  the  (100) 
dislocations  in  CoZr  (as 
well  as  in  CoTi)  have  a 
much  wider  core  than  in 
NiAl  (Table  6).  Note,  that  the  observation  of  an  unusually  wide  core  of  (100)  dislocations  in 
CoTi  was  reported  previously  [39].  This  feature  of  the  dislocation  structure  will  result  in  a 
low  Peierls  stress,  high  dislocation  mobility  and  hence  easy  plastic  relaxation  which  will 
prevent  brittle  fracture. 


0  0.1  0.2  03  0.4  0.5  0  0.1  0.2  0.3  0.4  0.5 

u/a  u/a 

Figure  14:  <100}{01 1 }  section  of  the  a-surface  in  (a)  CoZr  and  (b)  NiAl  before 
and  after  atomic  relaxation  perpendicular  to  the  slip  plane  direction. 


Thus,  the  high  ductility  of 
CoZr  is  an  intrinsic 
property  and  is  connected 
with  an  anomalously  low 
relaxed  GSF  energy  near 
1/2(100>{01 1 }  shear, 
which  is  also  an 
indication  of  the  phase 
instability  of  B2  CoZr.  A 

Figure  15:  Core  structure  of  the  <  1 00>{0 1 1 }  screw  dislocations  in  (a)  CoZr  detailed  understanding  of 
and  (b)  NiAl.  the  mechanism  of  this 

instability  will  require 

further  investigations.  Due  to  the  similarity  of  the  electronic  structure,  we  expect  similar 
behavior  for  all  Co-based  intermetallics. 
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10.  Dislocation  core  structure  in  B2  intermetallics  CoZr  and  NiAl:  large  scale  ARM  with 
ah  initio  parametrization. 

The  calculations  within  the  planar  PN  model  presented  above  allowed  us  to  identify  the 
important  features  of  the  dislocation  core  structure  in  CoZr.  However,  the  PN  model  could 
not  explain  the  temperature  yield  stress  anomaly,  which  is  related  to  the  existence  of  at  least 
two  (mobile  and  sessile)  core  structures;  therefore,  a  consideration  of  the  dislocation 
structure  beyond  the  scope  of  2D  PN  model  is  necessary. 
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To  make  clear  the  conditions  for  the  appearance  of  the  non-planar  dislocation  core  we 
investigated  the  structure  of  <100)  dislocations  in  B2  CoZr  and  NiAl  in  the  framework  of  a 
combined  approach  including  the  atomic  row  model  (ARM)  for  screw  dislocations  with  ab 
initio  parameterization  of  the  interatomic  row  interaction  energy.  This  approach  was 
previously  successfully  applied  to  modeling  the  alloying  effect  on  the  dislocation  structure  in 
bcc  Mo  and  explaining  the  main  trends  in  SSS  phenomena  [23],[29].  To  describe  the 
dislocation  structure  in  B2  intermetallics,  we  have  generalized  this  approach  by  taking  into 
account  (i)  the  geometry  of  the  packing  of  atomic  rows  in  B2  type  alloys  and  (ii)  three 
contributions  <J>Aa(u),  and  <J>ge(w)  to  the  interatomic  row  potential  originating  from 

the  interaction  between  atomic  rows  with  atoms  of  different  type  (A-B,  A-A  and  B-B, 
respectively). 


0  4  0  6 

u/a  u/a 

Figure  16:  (a)  The  atomic  row  displacement  (ARD)  energies  forNi.  Al-rows  (curves 
1,  2)  and  for  Co,  Zr-rows  (curves  1  ,  2  )  in  NiAl  and  CoZr;  (b)  unrelaxed  GSF  energy 
for  (100){001 }  and  (1 00){01 1 }  shear  (curves  1,  1  and  2,  2  ,  correspondingly)  in 
NiAl  (1.2)  and  CoZr  (1  ,2). 


16  for  CoZr  and  NiAl.  While  ARD  energies  for  Ni  and  Al  rows 
energy  for  Zr  row  is  much  smaller  in  comparison  with  the  Co 
similar  trend  takes  place  for  the  GSF  energy  in  {001 }  and  {011} 


The  dependencies  of 
®ab(u),  <S>aa(u),  and  <&bb(u) 
were  determined  by  fitting 
to  the  atomic  row 
displacement  energies 
and  AEl  along 
(001)  with  the  GSF  energy 
for  <100>{01 1 }  shear. 
Results  of  calculations  of 
A Ea  ,  A Eb  and  GSF 

energies  by  FLAPW  [4] 
method  are  shown  in  Fig. 
in  NiAl  are  close,  the  ARD 
row  in  CoZr  (Fig.  16(a));  a 
planes  (Fig.  1 6(b)). 
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Using  the  ARM  with  ab 
initio  determined  <J>Ab(u), 
^>aa(u),  and  <1 >BB(u) 
potentials,  we  found  that 
screw  dislocations  in 
NiAl  have  a  compact 
core,  which  agrees  with 
our  previous  theoretical 
results  [35]  and 
experimental  observation 
[41],  For  CoZr,  we 


Figure  17:  The  planar  (a)  and  non-planar  (b)  cores  of  the  (100)  screw  0^ta'nec^  an  unusual 
dislocation  in  CoZr.  result:  two  lyPes  of 

dislocation  structures 

with  planar  and  non-planar  cores  are  possible;  these  are  shown  in  Fig.  17  in  the  form  of 
relative  displacement  maps.  The  planar  core  structure  is  in  agreement  with  previous  ab  initio 
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PN  calculations  for  unrelaxed  GSF  and  correspond  to  the  splitting  of  the  (100)  dislocation 
according  to  the  b  =  b/2  +  b/2  scheme  on  the  {01 1}  plane.  A  non-planar  core  is  formed  by 
splitting  on  the  two  crossing  {011}  planes  which  results  in  a  loss  of  mobility  for  this  type  of 
dislocation.  The  non-planar  core  originates  from  highly  different  GSF  energies  in  {001 }  and 
{011}  planes  -  a  common  feature  for  all  Co-based  B2  intermetallics. 

Finally,  an  increase  in  temperature  can  lead  to  a  loss  of  dislocation  mobility  due  to  planar- 
non-planar  core  transformation  which  will  result  in  a  positive  oy(T)  dependence.  Thus,  our 
calculations  explain  this  feature  of  mechanical  behavior  of  B2  CoZr  and  allow  us  to  conclude 
that  the  yield  stress  anomaly  is  an  intrinsic  property  of  the  Co-based  intermetallic  alloys, 
rather  than  an  extrinsic  effect  as  was  previously  proposed  [42]. 
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